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Abstract 

A brief review of cosmology in some generalized modified gravity theories with vacuum polarization 
is presented. Stability question of de Sitter solution is investigated. 



DGP brane dynamics with vacuum polarization 

0^ , The acceleration universe force us to explain this phenomena and this stimulate an interest to a number 



' theories of modified gravity. One of such theories is the well-known Dvali-Gabadadze-Porrati (DGP) 

, braneworld model [HE], for instance, can lead to an accelerating universe without the presence of either 

a cosmological constant or some other form of dark energy. Generalizations of the DGP model can 
result in a phantom-like acceleration of the universe at late times, which is not excluded by observational 
data. We consider the simplest generic braneworld model with action of the form 



(1) 

<sJ '< Here, is the scalar curvature of the metric Gat hi the five-dimensional bulk, and R4 is the scalar 

^ I curvature of the induced metric gab — Gat — naUb on the brane, where is the vector field of the inner 

unit normal to the brane, which is assumed to be a boundary of the bulk space, and the notation and 
conventions of [3] are used. The quantity K = g°'^Kab is the trace of the symmetric tensor of extrinsic 
curvature Kab = g'^oSI c^b of the brane. The symbol L{gab, 4>) denotes the Lagrangian density of the 
^ I four-dimensional matter fields <j) whose dynamics is restricted to the brane so that they interact only with 

I the induced metric gab- AH integrations over the bulk and brane are taken with the corresponding natural 

volume elements. The symbols M and m denote the five-dimensional and four-dimensional Planck masses, 
respectively, A is the bulk cosmological constant, and A is the brane tension (it may be interpreted as 
cosmological constant on the brane). Note also that original DGP model apply A = 0. 

In Z2 symmetry case cosmological constant is equal from two side of the brane Ai = A2 = A, and 
Qv^ ' corresponding dynamical equation in FRW background looks like [H [S] : 
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^ , Here I — 2m? /M^^ H = a/ a is the Hubble parameter, and p is the matter energy density on the 

brane. Here and below, the overdot derivative is taken with respect to the cosmological time t on the 
brane. The expression under square root tend to zero during the evolution. It mean that in some time 
all time derivative from Habble parameter H , H, ... tends to infinity while H is finite. It's a new specific 
type of singularity and it was studied very well in the literature [6 . The term containing the constant 
C describes the so-called "dark radiation." We don't take into account this term, also we consider a 
spatially flat universe (fc = 0). The "±" signs in the solution correspond to two branches defined by the 
two possible ways of bounding the Schwarzschild-(anti)-de Sitter bulk space by the brane [71[S]. 

Now let us briefly consider dynamic of the DGP brane. Classical dynamics depends significantly on 
the bulk cosmological constant A and the full picture is sufficiently complicate so we refer you to the 
original paper focusing on only essential point for further statement (A > 0). The case A > is shown 
in FigHJ The graph of ((2]) in the (iJ^,ptot) plane in Fig [T] illustrates that in an expanding universe 
the matter density p decreases (except for a "phantom matter" which we do not consider in the present 
paper), and the point in the plane {H^, ptot) moves from right to left in Fig[TJ 
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Figure 1: Plot of relation ^ in the case A > 0. Case (a) corresponds to pc > 0, and case (b) corresponds 
to Pc < 0, where pc = 3m^ (A/6 — l^^)- 

A striking feature of Fig. [T] is that the value of the Hubble parameter in the braneworld can never 
drop to zero. In other words, the Friedmann asymptote — > is absent in our case. The upper and 
lower branches in Fig. [T] describe the two complementary braneworld models: branches AB and DB are 
associated with Brane 2 and Brane 1 of [9] , respectively, while branches AC and DC correspond to the 
lower and upper signs in Q , respectively, and describe the two branches with different embedding in the 
bulk. It should be noted that, in many important cases, the behaviour of the braneworld does not have 
any parallel in conventional Friedmannian dynamics (by this we mean standard GR in a FRW universe). 
For instance, the BC part of the evolutionary track corresponds to "phantom-like" cosmology with H > 0, 
even though matter on the brane never violates the weak energy condition. 

Now let us consider a some quantum effects. In general, quantum effects in curved space-time can 
arise on account of the vacuum polarization as well as particle production. It is well known that the latter 
is absent for conformally invariant fields (which we shall consider here) and that, in this case, quantum 
corrections to the equations of motion are fully described by the renormalized vacuum energy-momentum 
tensor which has the form [TUl 



(T;fe) = (m2/2880^2)(^.i^^^ _ 3^^^^ _ iy^fc^im^'" + J9^kR^) 

(3) 



+ (m3/28807r2)i(2i?.,;fe - 2g,kR''a ~ 2RR,k + y^kR^) 



where mi, m2 depend upon the spin weights of the different fields contributing to the vacuum polarization. 
This effect is known for a long time in cosmology. For example, it was demonstrate the possibility of 
singularity problem solution by some specific changes of m2 and ms . 

Since we work in flat FRW space-time, it is comfortable to use next relation: 

Pg^k2H'^ + k3{2HH + 6HH^ -H^). (4) 

Where parameters ^2 and take the next form [T^ [THl [H] : 

m2 _N + lliVi/2 + 62iVi + 1411iV2 - 28Nhd 



60(47r)2 60(47r)2 

m3 _ N + 6iVi/2 + 12iVi + 6nN2 - SNhd 
60(47r)2 ^ 60(47r)2 ' 



(5) 
(6) 



Here Ni is number of the fields with spin i contributing to the vacuum polarization: iV - a number of 
scalar fields, -/V1/2 - a number of fermion fields, A^i - a number of vector fields, N2 (= or 1) - a number of 
gravitons, Nhd - a number of conformal scalar fields. Note also that there is strong restriction on /c2 and 
fca in usual 4D-space-time. At least need fcs < with any k2 for de Sitter (and in particular Minkovsky!) 
solution stability. Since we see that Minkovsky solution is stable (there isn't particle production in 
conformal flat vacuum) we may use condition k^ < as observational data. 
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In order to assess the effects of tfie vacuum polarization on the dynamics of the braneworld, one must 
add pq to the matter density in ^ so that p — > p + in those equations. An important consequence of 
this operation is that the form of the equation of motion changes dramatically — the original algebraic 
equation changes to a differential equation ! The dynamical equation ([2|) now takes the form 



1 



HH=-H^ -3HH^ + i2k3) V-fca H'^ + im^H^ - ptot ± 3M " A/6 . (7) 



The goal of the present paper is to study the stability of the classical solutions when vacuum polar- 
ization terms are taken into account. The fc2-term in ([4]), which does not contain time derivatives of 
can only change the position of the future stable points. On the contrary, due to the fca-term in (jl]), 
some classical solutions can lose stability. Therefore, for simplicity (and without loss of generality), we 
set fc2 = in our calculations. If the brane has nonzero tension A, the stationary points of ([7]) in the 
case ^2 = can be found by substituting A into ^ and setting p = Q. After that, we linearize Eq. ^ at 
these stationary points and find the eigenvalues of the corresponding linearized system. The condition of 
stability of the stationary point is that its eigenvalues's real parts are negative. 

The eigenvalues at the stationary points where H — are given by 




= l[h± \//l' + 4/2j , (8) 

where we have made the notation 

h = -iH, (9) 

(10) 

Two different signs in Eq. (jlOp correspond to two different equations of motion, while, in Eq. ([S]), we 
have two different eigenvalues of a single equation. 

Since /i is negative, the eigenvalue yLt2 corresponding to the "— " sign in Eq. ([8|) is also always negative. 
Moreover, p,i is positive if and only if /2 is positive. As a result, the stability of a fixed point is equivalent 
to the condition /2 < 0. More careful investigation of this condition (for details see original paper [T5] ) 
show that part BC on the Fig. [T] corresponding to effective phantom behavior is unstable with respect 
vacuum polarization. 

General brane dynamics with vacuum polarization 

Now let us generalize obtaining result. We investigate any theory, which lead to dynamical equation in 
the next form: 

P^F{H), (11) 

where F{H) may be any algebraical function, which is don't contain time derivative of H . First of all 
we investigate the case fc2 = 0. Substituting in ()lip expression ^ for pq we may rewrite equation as 
dynamical system: 

^"^l (12) 
C=§^-3CH+^F{H)^f{H,C). ^ ^ 

Linearizing this system at the fixed point {Hq , 0) 

C-i%)oC+i§{^)oH, ^^^^ 
we find its eigenvalues 



'^^^^ = -2^(9b^o±^ig^)l+^g^U (14) 
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Since {§fj)o = —'^Hq < the real part of eigenvalues ^1^2 takes negative values if and only if 



dc 

df\ 1 fdF 



dHj„ 2Hoh\dHj^^^- ^^^^ 

So we can see that it need condition (f§)o > for stability. From another hand using (fTTj) we find 
(J^)q = ^ > or in this case the equivalent condition ^ > 0. By another words all regimes with 
^ < is unstable with respect to vacuum polarization. 

Rewriting last condition in the form q'^^q^ < we find that any regime with /? < (that is natural 

in expanding universe) and with H > Q is unstable. Exactly similar regimes is called fantom-likc. The 
case describing here correspond for example to BC part in Fig. [TJ 

Now let us account the fc2 term contribute. It mean that a new term k2H^ is appear in equation 
pip . Transferring this term into right hand part of equation we may introduce a new function F'{H) = 
F{H) — k2H^ and all previous result is true for this function. So term k2H'^ is not infiuencing on the 
stability, but it may change dynamical equation and fantom regimes is disappear at all |16j . 



/(i?)-theories with vacuum polarization 

Now let us consider /(i?)-theories which more popular in resent time. The most general form of such 
theories may be written in the next form 17J (for a more general review of f{R) gravity, see aslo [T81IT9]): 

S^Srn + ^j d^xV^fiR), (16) 

here x = SttG and for the sake of simplicity we set 2% = 1 . The general equation of motion corresponding 
to HH]) is given by 

9tkf'f - f'n-k + f'Rik - ^fgik = Tik. (17) 
The trace of equation of motion p?)) reads 

3/;,' + f'R - 2/ = T. (18) 
So existence condition of de Sitter solution in vacuum is [501 HI] 

2/(i?o) - Rof'iRo) = 0, (19) 
and condition of its stability following from psp is given by 



RoriRo) 



- 1 > 0. (20) 



Now let us take into account the vacuum polarization effects. For the sake of simplicity we investigate 



only k2 — case. So we find from ^ (T/) = — Substituting this relation into (fT5)) we find that this 



equation is restored to a normal state (without (T/) contribute) by the / ^ f — j^R^ transformation. 
And substituting this one into (|20p we find new de Sitter stability condition: 

fiRo) - Rof"{Ro) 

i?o(/"(i?o) - f ) ' ^ ^ 

which is turn into the expression l|20p by the limit k^ = 0. For instance, let us consider meaning of new 
expression (^1]) for the model [52] f{R) = R + -R~™ + R" . We see that only denominator of the ((5T|) has 
change when vacuum polarization take into account, while numerator is unchangeable. From another 
hand for this class of theories /" — m{m + 1)R~™~^ + n{n + l)i?"~^ > and since ^3 < taking into 
account of vacuum polarization effect can't change denominator's sign. It mean that in this case quantum 
effects do not infiuence on de Sitter stability condition. Note also that vacuum polarization effects to 
above f{R) gravity model may suppress instabilities as it was noted in [23] . 
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Conclusion 



Some modified gravity theories was studied with respect to vacuum polarization. It was demonstrate 
instabiUty of effective fantom regimes in the brane cosmology caused by vacuum polarization. Modified 
condition of the vacuum de Sitter solution stability in f{R) theories has been derived. 
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